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A new transport equation for the intermittency factor was proposed to predict separated 
and transitional boundary layers under low-pressure turbine airfoil conditions. The in- 
termittent behavior of the transitional flows is taken into account and incorporated into 
computations by modifying the eddy viscosity, p t , with the intermittency factor, y. Tur- 
bulent quantities are predicted by using Menter’s two-equation turbulence model (SST). 
The intermittency factor is obtained from a transport equation model, which not only can 
reproduce the experimentally observed streamwise variation of the intermittency in the 
transition zone, but also can provide a realistic cross-stream variation of the intermittency 
profile. In this paper, the intermittency model is used to predict a recent separated and 
transitional boundary layer experiment under low pressure turbine airfoil conditions. The 
experiment provides detailed measurements of velocity, turbulent kinetic energy and in- 
termittency profiles for a number of Reynolds numbers and freestream turbulent intensity 
conditions and is suitable for validation purposes. Detailed comparisons of computational 
results with experimental data are presented and good agreements between the experi- 
ments and predictions are obtained. [DOT: 10.1115/1.1580159] 


1 Introduction 

Flow transition plays an important role in turbomachinery ap- 
plications. Majority of boundary layer flows in turbomachines in- 
volve flow transition under the effects of freestream turbulence, 
diverse pressure gradients, wide range of Reynolds numbers, flow 
separation, and unsteady wake-boundary layer interactions. 

Prediction of this type of complex flows is an important ele- 
ment in analysis and performance evaluation of gas turbine engine 
components and ultimately in the design of more efficient jet en- 
gines. Especially, in low-pressure turbine applications prediction 
of transition behavior is even more important for reasons of de- 
sign efficiency. For low-pressure turbines the flow is mostly tur- 
bulent at the high Reynolds number conditions encountered at 
take-off and the efficiency is at its design maximum. However, 
due to decrease of Reynolds number caused by high altitude con- 
ditions at cruise speeds, design based on the sea level conditions 
tends to underpredict losses and thus leads to substantial drops in 
efficiency [1—3]. These losses are attributed to flow separation on 
the suction surface of the turbine blades. At low Reynolds num- 
bers with low freestream turbulence, the boundary layers on the 
airfoil surface have a tendency to remain laminar and hence the 
flow may separate before it becomes turbulent. This may cause a 
drop in efficiency and result in increase of fuel consumption. The 
impact of such losses is directly felt on the operation costs. It has 
been estimated that a 1% improvement in the efficiency of a low 
pressure turbine would result in a saving of $52,000 per year on a 
typical airliner [4], 

In order to calculate the losses and heat transfer on various 
components of gas turbine engines, and to be able to improve 
component efficiencies and reduce losses through better designs, 
accurate prediction of development of transitional boundary layers 
is essential [1], For an accurate prediction of transitional flows 
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under the diverse conditions encountered in turbomachineiy ap- 
plications, the aim is to find a model providing the following 
properties: 

• Physically accurate and versatile: The model should be ca- 
pable of accurately predicting transitional flows under the diverse 
conditions encountered in turbomachinery applications, such as 
pressure gradients, freestream turbulence, wide range of Reynolds 
numbers, unsteady wake-boundary layer interactions and flow 
separation. 

• Computationally efficient and inexpensive: The model should 
not involve extensive computational effort and should be rela- 
tively cheap to compute. 

• Compatible with current CFD methods: The model should be 
easy to implement into existing CFD codes without requiring ex- 
tensive changes in computational strategy. 

One of the current methods for predicting flow transition is to 
use the stability theory. In this method, stability equations are 
solved at streamwise stations in order to predict the onset of tran- 
sition. This method requires prior solution of the mean flow field 
and returns only the onset point of transition without any infor- 
mation on the turbulent part of the flow. Another method is using 
empirical correlations in the form of e n . This type of methods also 
require prior solution of the mean flow field. These two methods 
are not compatible with current CFD methods. 

One method compatible with current CFD methods is the use of 
low- Reynolds number turbulence models. Savill [5,6], has orga- 
nized a number of workshops to assess the capability of current 
turbulence models in predicting flow transition. The comparisons 
showed that none of the existing models was adequate to predict 
flow transition for a range of flow conditions. Westin and Henkes 
[7] have also tested a large variety of turbulence models and com- 
pared models’ performances in predicting a few of the T3-series 
transition flow experiments [5]. They indicated that none of the 
models could predict both onset location and length of transition 
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for a range of flow conditions. This outcome is not at all surpris- 
ing since most of the current turbulence models are not designed 
to predict flow transition. 

An alternative method to this approach is to use the concept of 
intermittency to blend the flow from the laminar to the turbulent 
regions. This approach, although highly empirical, has shown 
some successes in predicting transition behavior. Dhawan and 
Narasimha [8] correlated the experimental data and proposed a 
generalized intermittency distribution function across flow transi- 
tion. The correlation was later improved by Gostelow et al. [9] for 
flows with pressure gradients subject to a range of freestream 
turbulence intensities. 

Solomon et al. [10], following the work of Chen and Thyson 
[11], developed an improved method to predict transitional flows 
involving changes in pressure gradients. In this model, the effects 
of changing streamwise pressure gradient on the breakdown phys- 
ics and spot spreading rates are taken into account. This is accom- 
plished by varying the spot spreading angle and propagation pa- 
rameter through the transition zone according to the local pressure 
gradient parameter. 

Steelant and Dick [12] proposed a transport equation for inter- 
mittency, in which the source term of the equation is developed 
such that the y distribution of Dhawan and Narasimha [8] across 
the transition region can be reproduced. Steelant and Dick used 
their model, coupled with two sets of conditioned Navier-Stokes 
equations, to predict transitional flows with zero, favorable, and 
adverse pressure gradients. However, since their technique in- 
volved the solution of two sets of strongly coupled equations, the 
method is not compatible with existing CFD codes, in which only 
one set of Navier-Stokes equations is involved. Moreover, the 
model was designed to provide a realistic streamwise y behavior 
but with no consideration of the variation of y in the cross-stream 
direction. 

Cho and Chung [13] developed a k-e-y turbulence model for 
free shear flows. Their turbulence model explicitly incorporates 
the intermittency effect into the conventional k-e model equa- 
tions by introducing an additional transport equation for y. They 
applied this model to compute a plane jet, a round jet, a plane far 
wake, and a plane mixing layer with good agreements. Although 
this method was not designed to reproduce flow transition, it pro- 
vided a realistic profile of y in the cross-stream direction. 

Suzen and Huang [14] improved the intermittency transport 
equation by combining the best properties of Steelant and Dick’s 
model and Cho and Chung’s model. Their model reproduces the 
streamwise intermittency distribution of Dhawan and Narasimha 
[8] and also produces a realistic variation of intermittency in the 
cross-stream direction. The model is capable of predicting flow 
transition under diverse conditions. The predicting capabilities of 
this model have been validated against T3-series experiments of 
Savill [5,6] and low-pressure turbine experiments of Simon et al. 
[15] with good success [14,16,17]. 

In the current research we focus on further validation of the 
intermittency transport model of Suzen and Huang [14] against 
the experiments of Hultgren and Volino [18], Hultgren and Voli- 
no’s experiments investigated the effects of freestream turbulence 
and Reynolds number on separated and transitional boundary lay- 
ers under low-pressure turbine airfoil conditions. In their experi- 
ments, a flat plate boundary layer subject to a streamwise pressure 
gradient was studied. The superimposed pressure gradients were 
produced by attaching a two-dimensional contoured shape to the 
wall opposite to the test surface and by applying suction on the 
contoured wall. The resultant pressure profile represents that on 
the suction side of the Pak-B airfoil. The experiments covered a 
range of flow conditions including Reynolds numbers between 
50,000 and 300,000 and freestream turbulence intensities between 
0.2 and 7%. These cases cover a realistic range of operating con- 
ditions from take-off to cruise. On the test wall, velocity, turbulent 
kinetic energy and intermittency profiles were measured at 14 
streamwise stations. In the measurements, quantities such as 


skin friction coefficients, transition start and end locations, and the 
locations of separation and reattachment were also determined. 
Further details of the measurements and experimental data were 
given by Volino and Hultgren [19]. 

The experiments of Hultgren and Volino [18] provide a good 
set of data for the development and validation of the models for 
flow transition. In the next section, the intermittency transport 
model is presented and implementation of the model is described 
along with the empirical correlations employed for the onset of 
transition. In Section 3, the numerical details of the prediction 
process are given. In Section 4, the predictions of the new inter- 
mittency model are compared against the experimental data. Fi- 
nally, conclusions are provided in Section 5. 


2 Transport Model for the Intermittency 

In this section, the transport model for intermittency is pre- 
sented. The model combines the transport equation models of 
Steelant and Dick [12] and Cho and Chung [13]. Details of the 
development and implementation of the transport model are given 
in Suzen and Huang [14,16], Suzen et al. [17]. 

The model equation is given by 
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The distributed breakdown function, /(.v) has the form 
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The shear stresses are defined as 
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The blending function F is constructed using a nondimensional 
parameter, k/Wv, where k is the turbulent kinetic energy and IF is 
the magnitude of the vorticity. The blending function has the form 


F= tanh 4 


k/Wv 


[ 200(1 — y 01 ) 0 - 3 

The model constants used in Eq. (1) are 


(5) 


a yt =(T yt = 1.0 C 0 =1.0 Ci = 1.6 

C 2 = 0.16 and C 3 = 0.15 

The intermittency is incorporated into the computations simply 
by multiplying the eddy viscosity obtained from a turbulence 
model, fx , , by the intermittency factor, y. Simon and Stephens 
[20] showed that by combining the two sets of conditioned 
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Navier-Stokes equations and making the assumption that the Rey- 
nolds stresses in the nonturbulent part are negligible, the intermit- 
tency can be incorporated into the computations by using the eddy 
viscosity, fi* , which is obtained by multiplying the eddy viscos- 
ity from a turbulence model, n , , with the intermittency factor, y; 
that is, 

AL* = 7AL (6) 

is used in the mean flow equations. It must be noted that y does 
not appear in the generation term of the turbulent kinetic energy 
equations. 

The intermittency model had been implemented into TURCOM 
code developed by Huang and Coakley [21] and validated against 
low-pressure turbine experiments of Simon et al. [15] in Suzen 
et al. [17], 

It is essential that the turbulence model selected to obtain fj. t 
must produce fully turbulent features before transition location in 
order to allow the intermittency to have full control of the transi- 
tional behavior. Menter’s [22] SST model satisfies this require- 
ment. It produces almost fully turbulent flow in the leading 
edge of the boundary layer and therefore it is used as a base- 
line model to compute ju, and other turbulent quantities in the 
computations [17]. 

The value of tier used in evaluating the constants given by Eq. 
(3) is provided by the following correlation for zero pressure gra- 
dient flows [17]: 

na=(nv 1 /U i )a= 1.8X 10 -n Tu 7/4 (7) 

When flows are subject to pressure gradients, the following 
correlation is used: 

no- | J i/(i-«P(o.75xio‘ic,Tu- < >’)) ) K (<0 

( n<r) ZPG ~ 10 — 322 < 5985 , £,>0 (8) 

with M defined as: 

Af=(850Tu~ 3 — 100 Tu~°' 5 + 120) 

where {n&) Z p C is the value for flow at zero pressure gradient and 
can be obtained from Eq. (7) and K,=(v/ U 1 t )(dUldx) t is the flow 
acceleration parameter. The favorable pressure gradient part of the 
foregoing correlation (for K t > 0) is from Steelant and Dick [12], 
The portion of the correlation for adverse pressure gradient flows, 
AT ( <0, is formulated using the transition data of Gostelow et al. 
[9] and Simon et al. [15] in Suzen et al. [17], 

The current approach uses the intermittency transport model to 
obtain the intermittency distribution for the transitional flows, 
while the onset of transition is defined by the correlations dis- 
cussed in the next section. 

Attached-Flow Transition. Abu-Ghannam and Shaw [23] 
suggested that the onset of transition for attached flows can be 
obtained by correlating the boundary layer momentum thickness 
Reynolds number to the freestream turbulence intensity according 
to 

F(\ e ) 

Re0t = 163 + exp F(\ e ) (9) 

where 

F(\ e ) = 6.91+12.75\ 9 +63.64X4 for \„<0, or 
F'(\ s ) = 6.91 + 2.48\ 9 — 12.27A.0 for \ s >0 

Although Abu-Ghannam and Shaw correlation shows good 
agreement with experimental data for flows with zero and adverse 
pressure gradients, the model is not very sensitive to flows subject 
to strong favorable pressure gradients, in which one would expect 
the transition to be delayed as a result of flow acceleration [17], 

To allow for a more sensitive response to strong favorable pres- 
sure gradients while maintaining the good features of Abu- 
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Fig. 1 Onset of separated flow transition 


Ghannam and Shaw in adverse pressure gradient region, the tran- 
sition criterion was re-correlated to the freestream turbulence 
intensity, Tu, and the acceleration parameter, K , , 

Re e ,= (120+ 150Tu“ 2/3 )coth[4(0.3— r,X 10 5 )] (10) 

where K, was chosen as the maximum absolute value of that 
parameter in the downstream deceleration region [17]. Equation 
(10) was designed to have a better fit of the available experimental 
data: while the correlation fits the transition data well for flows 
under adverse pressure gradients, it was purposely designed to rise 
rapidly as K, becomes positive. This measure is to reflect the fact 
that the flow becomes less likely to have transition when subject 
to favorable pressure gradients. 

Separated-Flow Transition. Roberts [24] proposed a semi- 
empirical theory to predict onset of transition within a laminar 
separation bubble over the airfoil suction surface. The transition 
Reynolds number, Re„ , which is based on the length defined be- 
tween the onset location of separation and that of transition, is 
correlated as a function of a turbulent factor representing effects 
of the external turbulence level and its disturbance spectrum. The 
model was simplified by Davis et al. [25] to only a function of 
turbulence intensity, 

Re J( =2.5X 10 4 log] o coth( 0. 1 732Tu) (11) 

where Tu is the freestream turbulence intensity value at the onset 
of separation. Although this model was originally proposed for 
swept wing flows, it had been used widely for predicting onset of 
transition in a variety of separated flows. However, comparisons 
of recently available data for onset of transition in separated flows 
suggest that the correlation for the onset of transition is better 
represented by a function of more than one parameter. 

In Fig. 1, Davis et al. [25] correlation is given along with the 
data of Roberts [24] (the original data used to develop Davis’ 
correlation), Simon et al. [15], Hultgren and Volino [18] and Lake 
[26]. Each point (except for the data of Roberts [24]) is numbered 
to indicate the Reynolds numbers of the cases based on exit ve- 
locity and suction surface length. There is some scatter in the data 
especially for high Reynolds number cases of Lake [26]. In order 
to capture the scatter of the data, we incorporated the effect of 
Reynolds number into the correlation for Re sr . The new correla- 
tion expresses Re s , in terms of the turbulence intensity (Tu) and 
the momentum thickness Reynolds number at the point of sepa- 
ration (Re&) and is given as; 

Re s( = 874 Re^ 71 exp[ - 0.4Tu] (12) 
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Fig. 2 Onset of separated flow transition with new correlation 
(Eq. 12) 


In Fig. 2, the new correlation is compared with the experimen- 
tal data for different values of Re ft . . In contrast to Davis et al. [25] 
correlation, the current formula seems to provide a better repre- 
sentation of the data scattering. In the present paper Eq. (12) is 
used as a replacement of Davis et al. correlation to predict the 
onset of separated-flow transition. 

3 Numerical Details 

As mentioned earlier, the experiments made use of a flow suc- 
tion on the upper (contoured) surface to arrive at the desired pres- 
sure distribution. Although the shape of the upper contoured sur- 
face is known, no information is available for the conditions of the 
suction on this surface. This somewhat makes the computation of 
the problem ill-posed. In order to provide the upper boundary of 
the computational domain, we defined a prescribed streamline, 
obtainable by integrating the experimental measured velocity pro- 
files, as the upper boundary. Once the streamline is defined as the 
upper computational boundary, a slip boundary condition can be 
applied at the upper surface. 

For example, in Fig. 3 the streamlines obtained by using differ- 
ent values of 



Fig. 3 Streamlines from experimental data, Re=300,000, 
FSTI=7% 
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are presented for the case with Re=300,000 and FST1=7%. In 
this figure 'Pq is the 'P obtained from Eq. (13) at the first station 
taking h as the maximum height of the measured y position. To 
define the streamlines, we performed the integration of measured 
velocity profiles based on Eq. (13) in all x stations such that 'P is 
a constant. Any one of the streamlines shown in Fig. 3 can be used 
as the outer boundary of the computational domain, as long as 
they are well outside the boundary layer. Similar streamline 
curves can be generated for other cases. 

In the computations, we selected a streamline corresponding to 
'P = O.V'Fo as the upper boundary of the computational domain. 
The shapes of the upper streamlines are presented in Fig. 4 for all 
of the cases considered. In this figure the symbols show the points 
obtained from integration of the experimental velocity profiles and 
the lines indicate the boundary shapes used in the computations. 
By varying different streamline positions, we have found the 
streamline corresponding to 'P = 0.7^0 is sufficiently remote 
from the boundary layer to be considered as freestream. The 
boundary layer edge at each axial station is taken to be the loca- 
tion where the velocity is 99% of the maximum streamwise ve- 
locity at that station. Nonuniform meshes consisting of 140X60 
grid points are used for all calculations. All grids have first y + 
values less than 0.5. A grid refinement study for the case with 
Re=50,000 and FSTI=7% (using a 250X100 mesh) has revealed 
that the solution obtained by the current mesh can be accepted as 
a grid independent solution. 

Although the last experimental station is located at x/L = 1 .06, 
the computational domain is extended to x/L= 1.4 in the stream- 
wise direction to allow the outflow boundary condition to be ap- 
plied in the computations. We have extended the upper streamline 
boundary up to x/L = 1.25 and a constant pressure outflow bound- 
ary conditions is applied downstream of that station. Although this 
extension of the domain may appear to be arbitrary, it is necessary 
to avoid having reversed flow close to the exit boundary. Since no 
reverse flow was encountered for all cases after x/L*= 1.25, the 
extension of the domain up to x/L=\A seems reasonable. The 
plate surface is assumed to be an adiabatic wall and a constant- 
pressure outflow boundary condition is applied on the exit plane 
for all cases. 

In order to obtain accurate inlet profiles for the computations, a 
laminar computation over a flat plate is performed. From this 
computation, the profile matching the momentum thickness the 
first experimental station is extracted as the inlet condition for the 
computations. 
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Fig. 5 Overall view of effects of Reynolds number and freestream turbulence intensity 


4 Results and Discussion 

In this section the comparison of the predictions with experi- 
mental data of Hultgren and Volino [1 8] are presented. 

The experiment considered in this study involves two different 
Reynolds numbers, Re=300,000 and Re=50,000 (based on nomi- 
nal exit velocity, U re f, and wetted plate length, L = 0.208 m) and 
two freestream turbulence intensities, FSTI=7% and FSTI=0.2%. 
Figure 5 depicts the dynamic interplay between transition and 
separation when subject to variations in Reynolds number and 
freestream turbulence effects. In this figure the insets s and t de- 
note the onset of separation and transition, respectively, and the 
inset r denotes the reattachment point. The first case considered is 
with Re=300,000 and FSTI=7%. For this case, the transition 
position is at x/L = 0.66 and only a tiny flow separation is ob- 
served near x/L*** 0.7. As the Reynolds number is decreased from 
300,000 to 50,000 while keeping the freestream turbulence inten- 
sity the same, the onset of transition is delayed until x/L = 0.85. 
As a result of the delay of the transition, a sizable flow separation 
is encountered between x/L*** 0.65 and .0. On the other hand, 
when the turbulence intensity is decreased from 7 to 0.2% while 
the Reynolds number is maintained at 300,000, the transition on- 
set is delayed until x/L*** 0.79. As a result, a small but visible 
separation bubble is observed in the vicinity of x/L*** 0.75. The 
last case involved a simultaneous decrease of Reynolds number 
and freestream turbulence intensity, from 300,000 to 50,000 and 
from 7 to 0.2%, respectively. In this case the onset of transition is 
delayed to x/L*** \. 03 and a massive separation extending from 
x/L ***0.65 to ^1.25 is observed. The current prediction is an at- 
tempt to mimic the aforementioned dynamic behavior of the in- 
teraction between transition and separation when subject to 
changes in Reynolds number and freestream turbulence intensity 
conditions. 

For the case involving Re=300,000 and FSTI=7%, experimen- 
tal onset point of transition is at Re s =336 which corresponds to a 
location x/L = 0.66. ITowever, the onset point of transition for this 
case seems to be higher than other experimental data trend com- 
piled by Mayle [1] and Savill [5] and also higher than the value 
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given by the correlation of Abu-Ghannam and Shaw [23]. For 
example, the freestream turbulence intensity at the point of tran- 
sition is Tu=1.7% and the acceleration parameter is, K,= - 2.5 
X10~ 6 , by applying Eq. (10), the predicted onset location is 
much earlier (at Re e =230, corresponding to x/L = 0.55). To show 
the effect of the transition positions, we have performed two dif- 
ferent computations for this case: one using the experimental on- 
set point of transition and the other using the onset point obtained 
from the correlation, Eq. (10). 

The predicted pressure coefficient distribution along the surface 
is compared with experimental data, as shown in Fig. 6(a). In this 
figure, the results of the computation using the experimental onset 
of transition point is titled as “Computation 1” and the one ob- 
tained by utilizing Eq. (10) is denoted as “Computation 2.” Both 
results show very good agreements with experimental data and 
only slight differences in the distributions near x/L*** 0.65 are ob- 
served between the two computations. 

The velocity profiles for the two computations are compared 
with experimental data in Figs. 6(b) and (c). In the experiment, the 
flow remains laminar before station x/L = 0.69. At station x/L 
= 0.69 the velocity profile indicates that the flow is on the verge 
of separation. At the next station (x/L = 0.75), the boundary layer 
is attached and the flow continues to develop as an attached tur- 
bulent boundary layer. Although it is not seen from the measured 
velocity profiles, it was reported experimentally that a very small 
separation region exists between the two measured stations, x/L 
= 0.69 and 0.75. In contrast, no separation was observed in both 
computations although the comparison of the velocity profiles 
shows the predictions agree very well with experimental data. 
Both computations produce similar results except in the region 
x/L ***0.7. Due to the fact that the onset point of transition of 
computation 2 is slightly upstream, the result of computation 2 is 
fuller and is less likely to separate than that of computation 1; 
otherwise, the differences between the two profiles are small. 

The intermittency profiles of the two computations are com- 
pared with experimental data in Figs. 7(a) and (ft). As can be seen 
from Fig. 7a, the flow is laminar at the first seven stations. In both 
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Fig. 6 Prediction of experiment of Hultgren and Volino [18]— 
Re=300,000, FSTI=7% 


the experiment and the computation 1 , transition begins at station 
x/L = 0.66 and in computation 2 the transition starts at x/L 
= 0.55. While the transition of the computation 2 is ahead of the 
computation 1 , the length of transition is somewhat longer due to 
the fact that a small value of spot generation rate, ncr, is calcu- 
lated. The maximum value of the intermittency in the profiles of 
computation 1 reaches unity near station x/L = 0.81 whereas in 
computation 2, the intermittency reaches unity only after x/L 
= 0.88. Overall, computation 1 shows slightly better agreement 
with experiments than computation 2 even though the differences 
of the two results are mainly limited to the near wall region. 

In Figs. 8(a) and ( b ), comparisons of the turbulent kinetic en- 
ergy profiles are made. As can be seen from these figures, the 
freestream decay of turbulent kinetic energy throughout the plate 
is made to match the decay of freestream turbulence. This match 
provides the initial conditions for dissipation of the turbulent ki- 
netic eneigy (see Suzen and Huang [14], for details). The pre- 
dicted profiles of the turbulent kinetic energy seem to follow the 
experimental trend well even though the values are slightly larger 
before x/L = 0.7 and smaller after that. Again, the differences of 
the two computations were only limited to region near x/L = 0.7. 

The next case considered involves the same freestream turbu- 
lence intensity (FSTI=7%) while Reynolds number is reduced 
from 300,000 to 50,000. In the experiment, the flow separated at 
x/L = 0.63 and the onset of transition was observed at x/L 
= 0.85. In the computations, separation took place at x/L = 0.66 
and the onset point of transition calculated from Eq. (12) was at 
x/L = 0.88. The pressure coefficient distribution is compared with 
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Fig. 7 Comparison of intermittency profiles along the 
surface— Re=300, 000, FSTI=7% 


experimental data, as shown in Fig. 9(a). The agreement with the 
experimental data is good except in the region after x/L = 0.95, 
where it can be seen that the pressure coefficient is slightly un- 
derpredicted indicating that the streamwise velocity distribution is 
slightly overpredicted. 

The comparisons of the velocity profiles are given in Figs. 9(b) 
and (c). The profiles are in good agreement with the experiment 
up to x/L = 0.88. The experimental data shows that the flow at- 
taches after x/L = 0.94 whereas the predictions indicate that the 
separation is extended up to x/L = 1 .0. 




Fig. 8 Comparison of turbulence intensity profiles along the 
surface— Re=300, 000, FSTI=7% 
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Fig. 9 Prediction of experiment of Hultgren and Volino [18]— 
Re=50,000, FSTI=7% 



Y 


Fig. 10 Comparison of intermittency profiles along the 
surface — Re=50,000, FSTI=7% 


can be seen that the prediction gives rise to larger values of the 
turbulence kinetic energy for flow up to x/L*** 0.85. 

Finally, both Reynolds number and the ffeestream turbulence 
intensity are decreased to 50,000 and to 0.2%, respectively. For 
this case, flow separation was observed at x/L = 0.63 and the onset 
of transition was observed between stations x/L = 1 .0 and 1 .06 in 
the experiment. The flow did not reattach at the last measured 
station (x/L = 1.06). The predictions indicate that the flow sepa- 
ration is at x/L = 0.64 and the onset of transition calculated from 
Eq. (12) is sAx/L= 1.03. The predicted pressure coefficient distri- 


The intermittency profiles are compared to experimental data in 
Figs. 10(a) and ( b ). The general profile trend is well predicted 
even though the exact profile shapes are not well captured. In 
Figs. 11(a) and ( b ), comparisons of turbulence intensity profiles 
are shown. It can be seen that the comparisons between the pre- 
dictions and experimental data are favorable. 

The next case involves FSTI=0.2% and Re=300,000. In the 
experiment, the flow separated at x/L = 0.67, and the transition 
onset was observed between stations x/L = 0.75 and x/i = 0.81. 
The prediction shows that separation was at x/L = 0.66 and the 
onset of transition (obtained from Eq. (12)) was atx/L = 0.79. The 
comparison of the pressure coefficient distribution is presented in 
Fig. 12(a). The computed profile shows very good agreement with 
the experimental data. 

The comparison of the velocity profiles are given in Figs. 12 (b) 
and (c). The agreement between the predictions and experimental 
data is very good. It should be noted that there seems to be some 
discrepancies between the prediction and the measurement in the 
near-wall region near the flow separation region (x/L**=0.$), this 
difference is caused by the failure of the hot wire measurement for 
the flow reversal. 

The intermittency profiles are compared to experimental data in 
Figs. 13(a) and (b). The trend for the streamwise development of 
the intermittency factors seems to be well predicted, even though 
the model predicts a less diffusive behavior of the intermittency 
profiles near the ffeestream region. 

The turbulence intensity profiles along the surface are com- 
pared to experimental data in Figs. 14(a) and (b). Although the 
general trend of the turbulent kinetic energy profiles is captured, it 



Fig. 11 Comparison of turbulence intensity profiles along the 
surface— Re=50, 000, FSTI=7% 
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Fig. 12 Prediction of experiment of Hultgren and Volino [18]— 
Re=300,000, FSTI=0.2% 


bution is compared with experimental data, shown in Fig. 15(a). 
The agreement between the prediction and experimental data is 
very good. 

The comparisons of the velocity profiles are shown in Figs. 
15(6) and (c). As can be seen from the figures, good agreement is 
also observed for the velocity profiles. Once again, the discrep- 
ancy between the prediction and experiments near the wall is 
caused by the failure of the hot wire measurement in the separa- 
tion region. 

The intermittency profiles are compared to experimental data in 
Figs. 16(a) and (6). As can be seen from the figures, the flow 
remains laminar up to station x/L = 1.0. At the last measurement 
station, x/L= 1.06, the predicted intermittency factors show a 
larger magnitude than the experiments. 

The turbulence intensity profiles are compared to experimental 
data in Figs. 17(a) and (6). It can be seen that the prediction gives 
rise to larger values of the turbulent kinetic energy profiles for 
region jc/L>0.65. Even though the agreement is not that good, 
due to the fact that the transition only occurs at x!L= 1.03 its 
impact to the velocity profiles is not that pronounced. 

5 Concluding Remarks 

A new transport equation for the intermittency factor is em- 
ployed to predict a recent transitional boundary layer flow experi- 
ment under low pressure turbine airfoil conditions. The intermit- 
tent behavior of the transitional flows is taken into account by 
modifying the eddy viscosity with the intermittency factor. The 
new transport model not only can reproduce the experimentally 
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Fig. 13 Comparison of intermittency profiles along the 
surface— Re=300, 000, FSTI=0.2% 


observed streamwise variation of the intermittency in the transi- 
tion zone, but also provides a realistic cross-stream variation of 
the intermittency profile. Computations are performed for two dif- 
ferent Reynolds numbers and two different values of free stream 
turbulence intensities. Detailed comparisons with experiments are 
made for pressure coefficients, velocity, intermittency and turbu- 
lent kinetic energy profiles. Overall, good agreement with the ex- 
perimental data is obtained. It has been demonstrated that the 
predictions accurately mimic the dynamic behavior of the inter- 



Fig. 14 Comparison of turbulence intensity profiles along the 
surface— Re=300, 000, FSTI=0.2% 
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Fig. 15 Prediction of experiment of Hultgren and Volino [18]— 
Re=50,000, FSTI=0.2% 



Y 


Fig. 16 Comparison of intermittency profiles along the 
surface— Re=50, 000, FSTI=0.2% 
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Fig. 17 Comparison of turbulence intensity profiles along the 
surface— Re=50, 000, FSTI=0.2% 


play between the transition and separation when subject to varia- 
tions of Reynolds number and freestream turbulence intensity 
conditions. 

Acknowledgments 

This work is supported by NASA Glenn Research Center under 
grant NCC3-590. The project is part of the Low Pressure Turbine 
Flow Physics program of NASA-Glenn. This paper was originally 
published as AIAA Paper 2001-0446. 

Nomenclature 

C p = pressure coefficient, \—(\J e tU re j) 2 
FSTI = freestream turbulence intensity (%) 

K, = flow acceleration parameter, (v/ U 2 )(dU/ds) 
k = turbulent kinetic energy 
L = nominal suction surface wetted length 
N = nondimensional spot breakdown rate parameter, 
ncrO^/v 

n = spot generation rate 
p = static pressure 
Re = Reynolds no., L L / ref / v 
R e J( = (s,-s s )UJv 
Re*, = 0,U e lv 

s = distance along streamwise coordinate 
Tu = local turbulence intensity (%), u' !U 
U = boundary layer streamwise velocity 
U e = local freestream velocity 
U Kf = nominal exit freestream velocity 
u T = friction velocity 
W = magnitude of vorticity 

x = wetted streamwise distance along suction surface 
y n = distance normal to wall 

y + = y n t* T iv 

y = intermittency factor 
0 = momentum thickness 
'Kg = pressure gradient parameter, (8 2 /v)(dU/ds) 
p = molecular viscosity 
p t = eddy viscosity 
v = pJp 
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V t = fl t /p 

p = density 

a = spot propagation parameter 

Subscripts 

e = fireestream 
s = onset of separation 
t = onset of transition 
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